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Abstract. In this work we present some applications about the use of the so-called 
Cosmography with GRBs. In particular, we try to calibrate the Amati relation by using 
the luminosity distance obtained from the cosmographic analysis. Thus, we analyze the 
possibility of use GRBs as possible estimators for the cosmological parameters, obtaining 
as preliminary results a good estimate of the cosmological density parameters, just by using 
a GRB data sample. 
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1. Introduction 

It is a matter of fact that Gamma Ray Bursts 
(GRBs) are the most powerful explosions in 
the Universe; this feature makes them as one 
of the most studied objects in high energy 
astrophysics. The flux observed from their 
emission and the measurement of the red- 
shift z from the observations of the after- 
glow |Costa et al.( 1997)| point out a very high 
value for the isotropic energy emitted in the 
burst, so that there are some GRBs observed 
at very high redshift. Up to date, the far- 
thest GRB has a spectroscopic redshift of ~ 
8.2 |Tanvir et al.(2009)| [Salvaterra et al.(2009)| 
These interesting features suggest a possible 
use of GRBs as distance indicators; unfortu- 
nately our knowledge on the mechanisms un- 



derlying the GRB emission is not completely 
understood, so that their use as standard can- 
dles is still not clear. However, there exist 
some correlations among the observed spec- 
troscopic and photometric properties of the 
GRBs, allowing us to put severe constraints on 
the GRBs distances. What we need is an in- 
dependent estimation of the isotropic energy 
Ei S0 emitted from a GRB. Indeed, by using the 
GRB's fluence S t s measured by a detector in 
a certain energy range, it becomes possible to 
determine the luminosity distance d\ as follows 
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where Sboio is the bolometric fiuence emit- 
ted, obtained from the Schaefer formula 
|Schaefer(2007)l 



^l^ E * dE 



>bol 
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In literature there are many correlation formu- 
la{J each of them takes in account different 
observed quantities, but in this work we as- 
sume the validity of the so-called Amati rela- 
tion [AmSTetaL(2002)] for different reasons 

1 it relates just the isotropic energy with the 

peak energy in the vF(v) spectrum without 
considering others, 

2 it involves time-independent quantities, 

overcoming the problem of the large instru- 
mental biases, 

3 all of the long GRBs satisfy the Amati rela- 

tion, while the same is not true for the other 
correlations. 

In addition, although the Amati relation 
suffers of some biases, as the detector 
dependence of the observed quantity con- 
sidered I Shahmoradi & Nemiroff(2009)[ 



|Butler et al.(2009)| it seems to be well ver- 
ified from observations |Amati et al.( 2009) 
However, one of the most relevant challenge 
is represented by the calibration of the Amati 
relation, because a low redshift sample of 
GRBs is, up to now, lacking; a similar sample 
should be necessary in order to allow us 
to calibrate the relation too as well as the 
Supernovae la (SNela) calibration procedure. 
Anyway, a first computation of the relation 
parameters has been performed by consid- 
ering the concordance model, namely the 
Lambda-Cold Dark Matter (ACDM), obtain- 
ing a model-dependent luminosity distance. 
Unfortunately this procedure leads naturally 
to the so-called circularity problem when we 
take into account a cosmological use of the 
GRBs with the Amati relation. A possible 
solution has been provided by the use of SNela 
|Perlmutter et al.(1999)[ In other words, one 
can wonder if it is possible to calibrate GRBs 



by adopting at low redshift the SNela sample. 
This proposal has been already developed in 
literature by |Liang et al.(2008)] On the other 
hand, recently it has been investigated an al- 
ternative to solve this controversy, by adopting 
a model-independent procedure described by 
Cosmography, which shall be clarified in the 
next section. 

2. The cosmographic Amati relation 

As stressed above, the necessity to ac- 
count a procedure which is based on a 
model-independent way for characterizing the 
Universe dynamics is essential; indeed, dif- 
ferent cosmological tests may be taken into 
account; unfortunately for any case, one of 
the major difficulty is related to choosing 
which may be considered the less model in- 
dependent one. One of these, first discussed 



by Weinberg Weinberg(1972) and recently by 
Visser Visser(2004) proposes to consider the 
waste amount of kinematical quantities as con- 
straints to discriminate if a model works well 
or not. Cosmography is exactly what we mean 
for that; we refer to it as the part of cosmol- 
ogy trying to infer the kinematical quantities 
as the expansion velocity, the deceleration pa- 
rameter and so on, just making the minimal 
assumption of a Friedman-Robertson- Walker 
(FRW) metrics, being ds 2 = c 2 dt 2 -a(t) 2 \dr 2 + 

r 2 (sin 2 OdO 2 + d<f> 2 ) ], |Weinberg(1972)} in par- 
ticular, it is based only on keeping the ge- 
ometry by assuming the Taylor expansion of 
the scale factor a{t). In this way we do not 
do predictions about the standard Hubble law, 
but only about its kinematical constraints; it is 
worth noting that once expanded as a Taylor se- 
ries the Hubble law it is consequent to expand 
the luminosity distance d\ too and then the dis- 
tance modulus pXz) Capozziello & Izzo(2010); 
unfortunately it is clear that a similar expan- 
sion diverges for z > 1. Thus to circumvent 
this mathematical issue it should be necessary 
to change the variable, defining conventionally 



y 



1 +z 



(3) 



For a review see Meszaros(2006) 



which limits the redshift range, i.e. y e (0,1). 
With this model-independent formulation we 
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can immediately determine the cosmographic 
parameters, in order to reconstruct the trend of 
the function di(y) also at high redshift. Indeed, 
our aim consists in assuming the luminosity 
distance obtained with a good distance indica- 
tors, (SNela), extending it also for high red- 
shifts. So far, as a first step we estimated the 
cosmographic parameters from a very large 
sample of SNela, by adopting the Union 2 
compilation Amanu flah et al.(2010))[ to per- 
form this, we used a likelihood function L cc 
e~ x t 2 , where the chi-squared, x 2 is given by 



X 



-z 



(4) 



Table 1. Cosmographic parameters obtained 
using the Snela sample Union 2. Note that we 
have considered for the determination of the 
jerk jo and of the snap sq the flatness condi- 
tion Clk = 0. The error on so does not include 
the contribute from covariance terms and the 
dimension of Hq are km (sMpc) . 



Parameter 


value 


error 


Ho 


69.90 


0.027 


qo 


-0.58 


0.03 


jo 


1.50 


0.22 


so 


-2.96 


1.58 



where yu, are the distance modulus for each 
Union SNela and cr/i, its correspondent error. 
The results are summarized in Table Q] 

Once having an expression for d/, in princi- 
ple, it would be possible to calibrate the Amati 
relation too, by using the observed redshift and 
the bolometric fluence Sbolo of a GRB, com- 
puting the isotropic energy, by inverting eq. 
[TJ Then, having as the Amati relation Ei S0 = 
A E y p ; , we evaluated the parameters A and y, 
through the use of a sample of 108 GRBs 



Capozziello & Izzo(2010) considering as esti 



mator a log-likelihood function and taking into 
account the possible existence of an extra vari- 
ability cr ext of the y data, due to some hid- 
den variables that we cannot observe directly 



D'Agostini(2005) The cosmographic calibra- 



tion gives as results the following values 

A = 49.17 ± 0.40, y= 1.46 ±0.29, (5) 

and in Fig.[T]the best fit curve in the E p - E, so 
plane is showed. 

3. Cosmological applications 

Although of its elegance, our calibration of the 
Amati relation has been obtained using a for- 
mulation of d/ which suffers from some theo- 
retical misleading problems. First of all, since 
it is defined for low values of the redshift, the 
consequent extension to higher redshift may 
bring to deviations from the real cosmological 
picture. In order to check this discrepancy, we 



plotted in fig. |2]both the distance modulus ob- 
tained from Cosmography by SNela and from 
a flat ACDM paradigm, with Q p = 0.27, be- 
ing the matter density. We immediately note a 
difference of one magnitude at redshift z»4, 
increasing with z, due to different possible rea- 
sons 

1 the propagation of the systematics in the 

analysis of the SNela used for calibration, 

2 the large scatter in the data sample of the 

Amati relation, 

3 the standard ACDM model fails at high red- 

shift. 

The latter assumption seems to be the less 
probable one, since the ACDM model is able 
to explain the growing of structure formations 
too. In addition, in the following, we are go- 
ing to present a cosmological application of the 
GRB sample to estimate the density parame- 
ters. 

Let us first compute the isotropic energy 
E, TO for each GRB from the cosmographic 
Amati relation, obtaining the distance modu- 
lus for each of them, by using the bolometric 
fluence Sboio of eq. ©. Thus, the GRB sam- 
ple becomes related to the following theoreti- 
cal distance modulus 



D L (z) = 



1+z 
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Fig. 2. Plot of the fj(y) computed for a fiducial 
ACDM cosmological model, the continuous line, 
and for the reconstructed fi(y) obtained from the 
cosmographic fit of the Snela, the dashed line, in 
function of the z redshift. 

with E(z) the reduced Hubble parameter, i.e. 
E = while Qjt represents the fractional 
curvature density at z — 0, and 

( sin(x), if £l k < 0, 
sinn(x) = < x, if Q.^ = 0, 

( sinh(x), if Q k > 0. 



We consider as likelihood L oc e x grbI 2 the 
function given by 

2 _ y 1 )t!h ~ /Jobs ' n\ 

where y. \, s is the observed distance modulus 
for each GRB, with error cry, ,-, derived from 
the Amati relation, while yu,/, is the value of 
the distance modulus evaluated from the cos- 
mological model. The constraints have been 
evaluated by a combined cosmological test, 
provided by the SNela, baryon acoustic os- 
cillations (BAO) and cosmic microwave back- 
ground (CMB). Hence the total x 1 is given by 
Wang & Mukherjee(2006)| 

X 2 = xIrb + Xsn + Xbao + X 2 cmb ■ (8) 

In order to perform it, we adopt the Union 2 
compilation Amanul lah et al.(2010))| deriving 
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the constraints and confidence limits by using 
the same statistic employed for GRBs. In par- 
ticular we adopt the CMB shift parameter R 



R = ^jn m Hlr(zcMB), 



(9) 



with 



r(ZcMB) 



in* 



-1/2 



sinn 



while the x 2 term is given by 
(R - Robs) 2 



XcMB ~ 



(10) 



where for R & s and its error we con- 
sider the recent WMAP 7-years observations 
|Komatsu et al.(2010)[ For the SDSS baryon 
acoustic oscillations (BAO) scale measurement 
and in particular the distance parameter A 



A = 
with 



r(ZBAo) 



2 czbao 



H{zbao) 



1/3 (Q m i/ 2 )^ 2 



CZBAO 



(11) 



r(ZBA0) = — 
"0 



- 1/2 sinn(mt| 1/2 



r 

Jo 



A BA0 = 0.469 (n,/0.98)- U j 5 and cr A = 
0.017 |Eisenstein et al.(2005)[ The redshift 
t<bao = 0.35 while the spectral index is 
reported in the respective technical paper 
|Komatsu et al.(2010)| The minimization of the 
total x 2 was done applying a grid-search 
method in the parameter space of the model 
considered. As a first analysis we considered 
again the case of the ACDM model, obtain- 
ing not good results, (see fig. |3). We conclude 
that this happened due to the lacking low- 
redshift GRB sample, so that we are not able 
to give a good accuracy for the best fit values 
obtained using the GRB data sample only. A 
natural extension of ACDM is represented by 
the wCDM model, the so-called Quintessence 
model Copeland et al.(2006) here again the 



results are not in good agreement with re- 
spect what we expected, (see fig. @). In or- 
der to show a good agreement with observa- 
tions we expect that, since GRBs are gener- 
ally at high redshift, a varying Quintessence 



Km 




Fig. 3. 68%, 95%, and 98% constraints on Q p and 
fi A in the ACDM model obtained from CMB (red), 
BAO (blue), the Union 2 Compilation (gray and 
blue) and the GRB sample considered in this pa- 
per (gray and black). The superimposed contour plot 
represents the combined final results. 



model, can provide the trend of the w- 
term, giving rise to a well-fitting procedure. 
Among all the possibilities we report below 
the so-called Chevallier-Polarski-Linder (CPL) 
IChevallier & Polarski(2001),Linder(2003)| as 



E(0 w(z) 



w + w a 



1+z 



(12) 



where w = j. In this way the distance modu- 
lus curve is sensitive to variations at high red- 
shift of the w quantity, and GRBs are the only 
source that can shed light on this topic. The 
performed analysis developed by using both 
the SNela and GRB data sample gives results 
quite in agreement with what we expected, see 
fig.[5]and tab. [2] (together with the other analy- 
sis) and seems to point out GRBs as fundamen- 
tal tracers of an evolving Dark Energy equation 
of state. 



4. Conclusions 

In this work we wondered if the possibility of 
using GRBs as distance indicators can be a real 
resource of the modern Precision Cosmology; 
obviously this deals with the issue that, up 
to now, we cannot claim that GRBs are stan- 
dard candles. We developed a statistical (com- 
bined) analysis in which the calibration of 
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Fig. 4. 68%, 95%, and 98% constraints on Q. p and 
w obtained from CMB (orange), BAO (green), the 
Union 2 Compilation (gray and blue) and the GRB 
sample considered in this paper (gray and black). 
The superimposed contour plot represents the com- 
bined final results. 
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Fig. 5. 68%, 95%, and 99.7% constraints on the 
CPL parameters w and w a obtained from the Union 
2 Compilation and the GRB sample. 



the luminosity distance has been performed 
by a SNela sample, testing different models 
(ACDM, wCDM and CPL parametrization) 
with a more complete sample, including GRB 
data. We obtain satisfactory results especially 
in the CPL case. We conclude that the present 
data cannot suggest to us something new about 
the standard model, but the procedure must be 
seen as a first application of the use of GRBs in 
cosmology, for future developments. In a next 
paper we shall present intriguing results, study- 



Table 2. Final results of our analysis on each 
cosmological model considered. 



Parameter 


value 


error 


ACDM 


"m 


0.279 


0.040 




0.726 


0.034 


n k 


-0.005 


0.001 


wCDM 




0.29 


0.08 


W 


-0.87 


0.15 


wCDM + CPL 




0.27 


0.08 


w 


-1.18 


0.38 




0.38 


0.59 



ing with more accuracy the quoted models and 
other alternatives. 
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